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Ballistic hot electron transport overcomes the well-known problems of conductivity and spin life- 
time mismatch that plagues spin injection in semiconductors with ferromagnetic ohmic contacts. 
Through the spin-dependent mean-free-path, it also provides a means for spin detection after trans- 
port. Experimental results using these techniques (consisting of spin precession and spin-valve 
measurements) with Silicon-based devices reveals the exceptionally long spin lifetime and high spin 
coherence induced by drift-dominated transport in the semiconductor. An appropriate quantita- 
tive model that accurately simulates the device characteristics for both undoped and doped spin 
transport channels is described; it can be used to determine the spin current velocity, diffusion con- 
stant, and spin lifetime, constituting a spin "Haynes-Shockley" experiment without time-of-flight 
techniques. A perspective on the future of these methods is offered as summary. 

KEYWORDS: spin polarized electrons, ballistic hot electron transport, spin injection and detection, 
semiconductor spintronics, spin precession 
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I. BACKGROUND 

The problems of spin-polarized electron injection and 
detection are central to the field of semiconductor spin- 
tronics. Since ferromagnetic metals have a large asymme- 
try between the spin-up and spin-down density of states 
and therefore high spin-polarization at the Fermi level, it 
is natural to assume that ohmic transport of these elec- 
trons into a semiconductor will provide a robust injection 
mechanism for spin-polarized currents. However, in the 
diffusive regime required by ohmic transport (where the 
electron mean free path (mfp) is smaller than all other 
lcngthscales and an electron temperature and chemical 
potential are well-defined), injection and detection effi- 
ciencies are suppressed by the large differences in con- 
ductivity and spin lifetimes between metals and semicon- 
ductors, resulting in the so-called "fundamental obstacle 
to spin injection" [J H i, [| . 

The problem can be simply viewed in the follow- 
ing way: Since the spin-up and spin-down conductivi- 
ties in the (non-magnetic) semiconductor are identical, 
a spin-polarized current must be driven by an asym- 
metric spin-dependent electrochemical potential (/if /J.) 
drop across the semiconductor. This necessarily requires 
Afi = /x-f — fj,^ =/= at the interface between the fer- 
romagnet (FM) and semiconductor. In terms of the 
semiconductor bulk conductivity as, the semiconductor 
spin transport lengthscale L, cross-sectional area A, spin- 
polarization P, and total current /, Ohm's law gives 



Afi = 2P ■ IR = 2PI 



(1) 



On the other hand, in the metallic FM the spin current 
flowing toward the FM/semiconductor interface is caused 
instead by a spin-dependent conductivity Py i- There- 
fore, the interfacial electrochemical splitting can also be 
written as 
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FIG. 1: Spin-dependent band diagram showing the electro- 
chemical potentials near the interface of a metallic FM (high 
conductivity and short spin lifetime) with a semiconductor 
(low conductivity and long spin lifetime). 



Afj, = I(Ri — -R-r) = / 
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where A is the lengthscale over which the non-equilibrium 
electrochemical splitting in the FM decays away from 
the interface. Since diffusive transport dominates in the 
metal because it cannot support internal electric fields, 
this lengthscale is the "spin diffusion length" s/Dt, where 
r is the spin lifetime in the FM metal and D is the dif- 
fusion coefficient. Unfortunately, this timescale is typi- 
cally very short, so the lengthscale A is also short. The 
spin-dependent electrochemical potentials near the FM- 
semiconductor interface are schematically illustrated in 

Fig. m 

With the definition of the FM bulk polarization as f3 = 
^ T ^ , we can rewrite the electrochemical splitting on the 
FM side (Eqn. [2]) as 
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where &fm is the average value of conductivity <TT ^ <T; . 
At an ideal interface where the absence of interfacial 
spin relaxationpreserves continuity of the electrochemi- 
cal potentials, 0, @| the expressions given for the interfa- 
cial electrochemical splitting in Eqns. [1] and [3] must be 
equal. This then allows us to calculate the polarization 
of the injected current in the semiconductor: 
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Note that the conductivity of FM metals is several or- 
ders of magnitude larger than any semiconductor, even 
in the degenerately doped regime. In addition, the spin 
diffusion length in the FM metal is many orders of magni- 
tude shorter than the spin transport length in the semi- 
conductor due primarily to the short (long) lifetime in 
the metal (semiconductor) . The polarization of the elec- 
tron current flowing through the semiconductor is then 
dependent on the product of two (very) small quantities, 
- 2 - s — and 4- ■ Unless the bulk conductivity of the FM is 
close to 1 (i.e. a 100% polarized ferromagnetic half-metal 
- still not possible in practice), ohmic injection of spin 
polarized electrons is therefore hopelessly inefficient. 

To circumvent this problem, the constraints of ohmic 
injection must be lifted. For instance, the requirement 
that the electrochemical splittings be equal on the FM 
and semiconductor side can be relaxed by inserting a re- 
sistive tunnel barrier, [ij The high conductivity of the FM 
can be addressed by using a ferromagnetic semiconductor 
which typically has a much lower carrier concentration 
than a metal. @, HI 

Here, we consider a mechanism for spin injection that 
differs fundamentally from ohmic injection in that 



1. Momentum scattering is not the smallest length- 
scale in the device; and 

2. Electrochemical potentials can not be uniquely de- 
fined in thermal equilibrium. 



The above-described obstacle to spin injection therefore 
does not apply. With our techniques, the physical length- 
scale of the metallic electron injection contacts is shorter 
than the mfp, and conduction occurs through states far 
above the Fermi level (as compared to the thermal en- 
ergy ksT), far out of thermal equilibrium. Because this 
transport mode utilizes electrons with high kinetic en- 
ergy that do not suffer inelastic scattering, it is known 
as "Ballistic Hot Electron Transport" . We will find that 
this transport mode is useful not only for spin-polarized 
electron injection into semiconductors, but also (and per- 
haps more importantly) for spin detection after transport 
through them. 
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FIG. 2: Schematic illustration of a tunnel junction used as a 
source for ballistic hot electron injection into a semiconductor 
conduction band. The emitter cathode electrostatic potential 
energy (qV) must exceed the Schottky barrier height q<f). 



II. HOT ELECTRON GENERATION AND 
COLLECTION 

Tunnel junctions (TJs, which consist of two metallic 
conductors separated by a thin insulator) [{J 

Ea nam 

0j 13 are the ideal hot electron source. The electro- 
static potential energy qV provided by voltage bias V 
tunes the energy of hot electrons emitted from the cath- 
ode, and the exponential energy dependence of quantum- 
mechanical tunneling assures a narrow distribution. Hot 
electrons thus created can be collected by the Schottky 
barrier if qV > q<fi. Because of its robust insulating 
native oxide, Aluminum (Al) is often employed during 
tunnel junction fabrication; although in principle any in- 
sulator can be used, it has been empirically found that 
the best are Al x O, MgO, and AlN.Qjl, [ll] The first pro- 
posal for a tunnel-junction hot electron injector was by 
Mead, who suggested another tunnel junction or the vac- 
uum level as a collector. [l7l [l8l [l9j Very soon thereafter, 
Spratt, Schwartz, and Kane showed that a semiconduc- 
tor collector could be used to realize this device with a 
Au/ALjO/Al tunnel junction. [20l| A schematic band dia- 
gram of this device is shown in Fig. O 

Metals have a very high density of electrons at and be- 
low the Fermi energy Ep. Therefore, sensitive hot elec- 
tron collection by the semiconductor conduction band 
relies on an electrically rectifying barrier to eliminate 
transport of these thermalized electrons across the metal- 
semiconductor interface, which would dilute the injected 
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hot electron current. This barrier is ideally created by 
the difference in work functions of the metal and the elec- 
tron affinity of the semiconductor [2ll I22I I23I HH, but in 
reality its height is determined more by the details of 
surface states which lie deep in the bandgap of the semi- 
conductor that pin the Fermi level. [25| There is of course 
always a leakage current due to thermionic emission over 
this "Schottky" barrier at nonzero temperature; because 
typical barrier heights are in the range 0.6-0.8 eV[26| for 
the common semiconductors Si and GaAs, hot electron 
collection with Schottky barriers is often performed at 
temperatures below ambient conditions to reduce current 
leakage to negligible levels. 



III. SPIN-POLARIZED HOT ELECTRON 
TRANSPORT 

The mfp of hot electrons in ferromagnetic metal films 
is spin dependent: majority ("spin up") electrons have 
a longer mfp than minority ("spin down") electrons. 
Therefore, an initially unpolarized hot electron current 
will become spin-polarized by spin-selective scattering 
during ballistic transport through a ferromagnetic metal 
film, where the polarization is given by 



P = 



(5) 



'•J +e A " 



where I is the FM film thickness, X ma j is the majority 
spin mfp, and \ m in is the minority spin mfp. This "bal- 
listic spin-filtering" effect can be used not only for spin 
polarization at injection, but also for spin analysis of a 
hot electron current at detection, much as an optical po- 
larizing filter can be used both for electric field polar- 
ization and analysis of photons by changing the relative 
orientation of the optical axis. The realization of this 
important spin analysis concept has enabled the spin in- 
jection, transport, and detection in Si described in this 
text. 



IV. SPINS IN SILICON 

Silicon is a relatively simple spin system|27]; Si could 
be considered as the "hydrogen atom" of semiconductor 
spintronics: 

1. A low atomic number (Z = 14) leads to a reduced 
spin-orbit interaction (which scales as Z 4 in atomic 
systems). 

2. As discussed previously, the preservation of spa- 
tial inversion symmetry of the diamond lattice leads 
to a spin-degenerate conduction band, eliminating 
D'yakonov-Perel' spin scattering in bulk Si. 



3. The most abundant isotope of Si (w 92%) is Si 28 , 
which has no nuclear spin, and this abundance can 
be fortified with isotopic purification to very high 
levels. [HI Therefore, the spin-nuclear (hyperfine) 
coupling is weak, as compared with other semicon- 
ductors where no such nuclear-spin-free isotope ex- 
ists. 

Due to its apparent advantages over other semiconduc- 
tors, many groups tried to demonstra te p henomena at- 
tributed to spin transport in Si[29l [30l l3ll l32l. l33l. 1341 l35l. 
HSH3 but this was typically done with ohmic FM-Si con- 
tacts and two-terminal magnetoresistance measurements 
or in transistor- type devices [38l l39j which are bound to 
fail due to the "fundamental obstacle" for ohmic spin in- 
jection mentioned in Sec. Hip], 0> B IU Although weak 
spin- valve effects are often presented, no evidence of spin 
precession is available so the signals measured are am- 
biguous at best. [13, Sl| Indeed, although magnetic ex- 
change coupling across ultra-thin tunneling layers of Si 
was seen, not even any spin- valve magnetoresistance was 
observed. [42 

These failures were addressed by pointing out that 
only in a narrow window of FM-Si Resistance-Area 
(RA) product was a large spin polarization and hence 
large magnetoresistance expected. [43[ Subsequently, sev- 
eral efforts to tune the FM-Si interface resistance were 
made. [3, [H, [H, H3, EU However, despite the ability to 
tune the RA product by over 8 orders of magnitude and 
even into the anticipated high-MR window (for instance, 
by using a low- work- function Gd layer), no evidence that 
this approach has been fruitful for Si can be found, and 
the theory has yet been confirmed only for the case of 
low-temperature-grown 5nm-thick GaAs.[49l. [50j 

Subsequent to the first demonstration of spin transport 
in Si,[51[ optical detection (circular electroluminescence 
analysis) was shown to indicate spin injection from a FM 
and transport through several tens of nm of Si, using first 
a Al^O tunnel barrier [52| and then tunneling through the 
Schottky barrier despite the indirect bandgap and 
relative lack of spin-orbit interaction which couples the 
photon angular momentum to the electron spin. These 
methods required a large perpendicular magnetic field 
to overcome the large in-plane shape anisotropy of the 
FM contact, but later a perpendicular anisotropic mag- 
netic multilayer was shown to allow spin injection into 
Si at zero external magnetic field. [54| While control sam- 
ples with non-magnetic injectors do show negligible spin 
polarization, again no evidence of spin precession was 
presented. 

More recently, four-terminal nonlocal measurements 
on Si devices have been made, for instance with AUO[55j 
or MgO[56j tunnel barriers, or Schottky contacts us- 
ing ferromagnetic silicide injector and detector. 57j Only 
Ref. l5fj| presents evidence of spin precession, but the low 
signal-to-noise obscures geometrical effects of precession 
expected from the oblique magnetic field configuration 
used. [Hi 
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FIG. 3: Schematic band diagram of a four-terminal (two for 
TJ injection and two for FM SMS detection) ballistic hot 
electron injection and detection device with a 350 /^m-thick 
Si transport layer. 



As mentioned in Sec. HI another way to overcome 
the conductivity/lifetime mismatch for spin injection 
is to use a carrier-mediated ferromagnetic semiconduc- 
tor heterointerface. The interfacial quality may have 
a strong effect on the injection efficiency so epitax- 
ial growth will be necessary. Materials such as dilute 
magnetic semiconductors Mn-doped Si[59(, Mn-doped 
chalcopyrites (6C" 



61, 62, 6£ 



64j , or the "pure" ferromag- 
netic semiconductor EuO [65[ have all been suggested, 
but none as yet have been demonstrated as spin injectors 
for Si. It should be noted that while their intrinsic com- 
pounds are indeed semiconductors, due to the carrier- 
mediated nature of the ferromagnetism, it is seen only in 
highly (i.e. degenerately) doped and essentially metallic 
samples in all instances. 



V. BALLISTIC HOT ELECTRON INJECTION 
AND DETECTION DEVICES 

Two types of tunnel junctions have been employed to 
inject spin-polarized ballistic hot electrons into Si. The 
first used ballistic spin filtering of initially unpolarized 
electrons from a nonmagnetic Al cathode by a ferro- 
magnetic anode base layer in direct contact to undoped, 
10-micron-thick single-crystal Si(100). [5l| Despite SVT 
measurements suggesting the possibility of 90% spin po- 
larization in the metal. [6a. \6l\ only 1% polarization was 
found after injection into and transport through the Si. 
It was discovered later that a nonmagnetic Cu interlayer 
spacer could be used to increase the polarization to ap- 
proximately 37%, (68| likely due to Si's tendency to read- 
ily form spin-scattering "magnetically-dead" alloys (sili- 
cides) at interfaces with ferromagnetic metals. 

Despite the possibilities for high spin polarization with 
these ballistic spin filtering injector designs, the short 
hot-electron mfps in FM thin film anodes causes a very 



small injected charge current on the order of lOOnA with 
an emitter electrostatic potential energy approximately 
leV above the Schottky barrier and a contact area ap- 
proximately 100 x 100/im 2 . Because the injected spin 
density and spin current are dependent on the product 
of spin polarization and charge current, this technique 
is not ideal for transport measurements. Therefore, an 
alternative injector utilizing a FM tunnel-junction cath- 
ode and nonmagnetic anode (which has a larger mfp) was 
used for approximately 10 times greater charge injection 
and hence larger spin signals, despite somewhat smaller 
potential spin polarization of approximately 15%.[69lr70j] 
These injectors can be thought of as one-half of a mag- 
netic tunnel junction. (7lT| with a spin polarization pro- 
portional to the Fermi-level density of states spin asym- 
metry, rather than exponentially dependent on the spin- 
asymmetric mfp as is the case with ballistic spin filtering 
described above. 

Although the injection is due to ballistic transport 
in the metallic contact, the mfp in the Si is only on 
the order of 10nm.[72| so the vast majority of the sub- 
sequent transport to the detector over a lengthscale of 
tens.JElt, WA , wdl l hundreds Wdi] or thousandsWk] of 
microns occurs at the conduction band edge following mo- 
mentum relaxation. Typically, relatively large acceler- 
ating voltages are used so that the dominant transport 
mode is carrier drift; the presence of rectifying Schot- 
tky barriers on either side of the transport region assures 
that the resulting electric field does nothing other than 
determine the drift velocit y o f spin polarized electrons 
and hence the transit time[74j, 1 79] - there are no spuri- 
ous (unpolarized) currents induced to flow. Furthermore, 
undoped Si transit layers are primarily used; otherwise 
band-bending would create a confining potential and in- 
crease the transit time, potentially leading to excessive 
depolarization (see Section IVl Bp . [761] 

The ballistic hot electron spin detector is comprised 
of a semiconductor-FM-semiconductor structure (both 
Schottky interfaces), fabricated using UHV metal- film 
wafer bonding (a spontaneous cohesion of ultra-clean 
metal film surfaces which occurs at room-temperature 
and nominal force in ultra-high vacuum) . [77| [78j After 
transport through the Si, spin-polarized electrons are 
ejected from the conduction band over the Schottky bar- 
rier and into hot electron states far above the Fermi 
energy. Again, because the mfp in FMs is larger for 
majority-spin (i.e. parallel to magnetization) hot elec- 
trons, the number of electrons coupling with conduction 
band states in a n-Si collector on the other side (which 
has a smaller Schottky barrier height due to contact with 
Cu) is dependent on the final spin polarization and the 
angle between spin and detector magnetization. Quan- 
titatively, we expect the contribution to our signal from 
each electron with spin angle 9 to be proportional to: 
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FIG. 4: In-plane magnetic field measurements show the "spin- 
valve" effect and can be used to calculate the spin polarization 
after transport. 
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Because the exponential terms are constants, this has the 
simple form oc cos 8+const; in the following, we disregard 
the constant term as it is spin- independent. 

The spin transport signal is thus the (reverse) current 
flowing across the n-Si collector Schottky interface. In 
essence, this device (whose band diagram is schematically 
illustrated in Fig. [3]) can be thought of as a split-base 
tunnel-emitter SVT with several hundred to thousands 
of microns of Si between the FM layers. 

Two types of measurements are typically made: "spin 
valve" in a magnetic field parallel to the plane of magneti- 
zation, and spin precession in a magnetic field perpendic- 
ular to the plane of magnetization. The former allows the 
measurement of the difference in signals between parallel 
(P) and antiparallel (AP) injector/detector magnetiza- 
tion and hence is a straightforward way of determining 
the conduction electron spin polarization, 



P = 



Ip — Lap 
Ip + Iap 



(7) 



Typical spin-valve measurement data, indicating «8% 
spin polarization after transport through 350/im undoped 
Si, is shown in Fig. [H 

Measurements in perpendicular magnetic fields reveal 
the average spin orientation after transit time t through 
the Si, due to precession at frequency uj — gfisB/h, 
where B is magnetic field. If the transit time is deter- 
mined only by drift, i.e. t — we expect our spin 
transport signal to behave oc cosg/isBt/h. However, 
due to transit time uncertainty At caused by random 
diffusion, there is likewise an uncertainty in spin preces- 
sion angle A9 — uAt which increases as the magnetic 
field (and hence ui) increases. When this uncertainty ap- 
proaches 2n radians, the spin signal is fully suppressed by 
a cancellation of contributions from antiparallel spins, a 
phenomenon called spin "dephasing" or Hanle effect [79j| . 
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FIG. 5: (a) A typical spin precession measurement shows the 
coherent oscillations due to drift and the suppression of sig- 
nal amplitude ( "dephasing") as the precession frequency rises. 
Our model simulates this behavior well, (b) The real part of 
the Fourier transform of the precession data in (a) reveals the 
spin current arrival distribution. 



On a historical note, our device is essentially a solid- 
state analogue of experiments performed in the 1950s 
that were used to determine the g-factor of the free elec- 
tron in vacuum using Mott scattering as spin polarizer 
and analyzer and spin precession in a solenoid [8(J. In 
our case, we already know the g-factor (from e.g. ESR 
lines), so our experiments in strong drift electric fields 
where spin dephasing is weak can be used to measure 
transit time with t = h/ ' g^isBi-K^ where Bi^ is the mag- 
netic field period of the observed precession oscillations, 
despite the fact that we make DC measurements, not 
time-of- flight [7-S. fT5l] . Typical spin-precession data, indi- 
cating transit time of approximately 12ns to cross 350/xm 
undoped Si in an electric field of » 580 V/cm, is shown 
in Fig. Oa). 

One important application of this transit time informa- 
tion from spin precession is to correlate it to the spin po- 
larization determined from spin- valve measurements and 
Eqn. [7] to extract spin lifetime. By varying the internal 
electric field, we change the drift velocity and hence aver- 
age transit time. A reduction of polarization is seen with 
an increase in average transit time (as in Fig. [Ha)) that 
we can fit well to first order using an exponential-decay 
model P tx e _t / T , and extract the timescale, r.[7Cj In this 
way, we have observed spin lifetimes of approximately 
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FIG. 6: (a) Fitting the normalized spin signal from in-plane 
spin-valve measurements to an exponential decay model us- 
ing transit times derived from spin precession measurements 
at variable internal electric field yields measurement of spin 
lifetimes in undoped bulk Si. (b) The experimental spin 
lifetime values obtained as a function of temperature are 
compared to Yafet's T~ 5,/2 power law for indirect-bandgap 
semiconductors [8l| and Cheng et al.'s T -3 derived from a 
full bandstructure theory. [82j 



VI. SPIN TRANSPORT MODELING 

Modeling the spin transport signal involves summing 
all the detector signal contributions of spins which be- 
gin at the injector at the same time and with the same 
spin orientation (i.e. initial conditions on spin density 
s(x, t = 0) = S(x)), where s(x, t) is determined by the 
drift-diffusion equation 

^£ = D— - v— - - (8) 
dt dx 2 dx r' 

where D is the diffusion coefficient and v is the drift 
velocity. The solution to this partial differential equation 
with the S(x) Dirac-delta initial condition is called the 
"Green's function", which can be used to construct the 
response to arbitrary spin injection conditions, including 
the DC currents used so far in experiments. 

Unlike open-circuit voltage spin detection [4ll l85l. l8rj||. 
our ballistic hot-electron mechanism employs current 
sensing. This must be accounted for in the model- 
ing, and since electrons crossing the metal-semiconductor 
boundary do not return via diffusion, it imposes an ab- 
sorbing boundary condition on the spin transport. [13] 
Whereas voltage detection is sensitive to the spin den- 
sity s(x = L,t) and the boundaries imposed on the drift- 
diffusion Green's function are at infinity, here our signal 
is sensitive to the spin current J s (x — L,t) = — D&\ x= l, 
where we must impose s(x = L,t) = 0. 

A. Undoped Si 



1/Lts at 60K in 350/xm-thick transport devices. 83] 1 . The 
temperature dependence of spin lifetime is compared to 
the T -5 / 2 power law predicted by Yafet [811] and the more 
recent (and more complete) theory of Cheng et al giving 
T- 3 [12 in Fig. Efb). 

If we can adequately model the transport and the ex- 
pected signals, we can make a spin-polarized electron ver- 
sion of the Haynes-Shockley experiment (originally used 
to measure diffusion coefficient, mobility, and lifetime of 
minority charge carriers ).[84j This experiment was very 
useful in the design of bipolar electronics devices such 
as junction transistors. By measuring the spin transport 
analogues of these parameters, we imagine enabling the 
design of useful semiconductor spintronic devices. 



1 In Ref. :70, a more conservative estimate of the spin lifetime (e.g. 
520ns at 60K) was obtained by fitting to the transit time depen- 
dence of an alternative quantity expected to be proportional to 
the spin polarization, rather than using Eqn. [7] directly. 



Assuming v is a constant independent of position x (as 
it is in undoped semiconductors), the Green's function of 
Eqn. [5] satisfying this boundary condition can be found 
straightforwardly using the method of images: 



s(x, t) = 



1 



2VnDt 



e mi 



, ,/n (a-2i,-«t)- 



The corresponding spin current at the detector is 



J s (x = L,t) = 



1 



L 



-t/r 



-t/r 



(9) 



(10) 



The Green's function in Eqn. [5] and the spin current 
derived from it in Eqn. [TU] (using typical values for v, D, 
and L at various times t) are shown in Fig. [7] 

Note that the expression in Eqn. [TU]is simply the s pin 
density in the absence of the boundary condition [85|, l86j 
multiplied by the spin velocity L/t. Because our mea- 
surements are under conditions of strong drift fields, we 
have approximated this term with a constant (e.g. the 
average drift velocity v) in previous work. [70j [8g, 189} To 
first order, this approximation involves a rigid shift of the 



7 




FIG. 7: Evolution of a Delta function spin distribution 
(Green's function) injected on the left side at time increments 
shown in the legend, with drift velocity 10 6 cm/s, diffusion co- 
efficient 100cm 2 /s, and absorbing boundary conditions at the 
detector (x = 100/im). Inset: Simulated spin current, given 
by the gradient of the Green's function at the detector, as in 
Eqn. [TO] 



signal peak to lower arrival times of ss 2D/v 2 and a sub- 
sequent error in the measured velocity of 2D/L. There 
is no first-order error in the dephasing. Because typical 
values for these variables are on the order of v > 10 6 
cm/s, L > 10~ 2 cm, and 10 2 < D < 10 3 cm 2 /s, the rela- 
tive error associated with this approximation (2D/Lv in 
both cases) is then just a few percent. 

We can use this Green's function in the kernel of an 
integral expression to model, for instance, the expected 
precession signal as the weighted sum of all the cosine 
contributions from electrons with different arrival times 
(0 < t < oo): 



J s (x — L,t) cos Lutdt . 



(11) 



FIG. 8: Depletion-approximation conduction-band diagrams 
of a 3.3pm-thick n-type (7.2 x 10 14 cm -3 ) doped Si spin- 
transport layer with injector-detector voltage drop of 3V 
(where the transport is dominated by diffusion against 
the electric field at the detector side and results in the 
exponential-like arrival distribution with typical timescale 
1000ns shown in the top inset); 4.46V (where the bias is 
enough to eliminate the electrostatically neutral region and 
fully deplete the transport layer; and 6V (where the potential 
well has been eliminated and transport is dominated by drift 
in the unipolar electric field, resulting in the gaussian-like ar- 
rival distribution with typical timescale of 0.050ns shown in 
the bottom inset). 
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Note that because J s is causal (i.e. J s {t < 0) = 0), we 
can extend the lower integration bound to —00 so that 
this expression is equivalent to the real part of the Fourier 
transform of J s (x — L,t). We can therefore use preces- 
sion measurements as a direct probe of the spin current 
transit time distribution by making use of the Fourier 
transform, as shown in Fig. [5jb). By fitting the data to 
the simulation, we can extract the transport variables v 
and D; along with the spin lifetime measurement, this 
comprises a spin Haynes-Shockley experiment. 

Eqn. [Tl] is valid only for the case of a purely perpen- 
dicular magnetic field B = B z z.^^ In the most general 
case, if there is also an in-plane field component B y y (due 
to misalignment or a static field) and/or if the perpen- 
dicular component of the field B z is strong enough to 
overcome the in-plane magnetic anisotropy of the film, 
the cosut term must be replaced with 



where in the simplest case the magnetization rotation 
angles #1,2 = arcsin(tanh (B z /rji,2)), and 771^ are the de- 
magnetization fields of injector and detector ferromag- 
netic layers (typically several Tesla).[76] From Eqn. [T21 
we can show that small misalignments from the out- 
of-plane direction during single-axis precession measure- 
ments cause only in-plane magnetization switching of in- 
jector and detector and a negligible correction to the am- 
plitude of precession oscillations to first order. [58( Fur- 
thermore, the expression is now composed of a coherent 
term proportional to cos lot and an incoherent term inde- 
pendent of t. In the case of extremely strong dephasing 
(caused e.g. by geometrically-induced transit length un- 
certainty), the first term contributes a negligible amount 
to the integral expression in Eqn. [11] and the signal is 
dominated by the field and magnetization geometry dic- 
tated in the second term.[73| 
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B. Doped Si 

In the above, we have assumed that the drift velocity 
v in Eqn. [5] was independent of x. This assumption is 
not valid in doped semiconductors due to the presence of 
inhomogeneous electric fields and band-bending caused 
by ionized dopants. Eqn. [5] then no longer has constant 
coefficients, so in general its Green's function cannot be 
solved analytically and a computational approach must 
be taken. 

Direct simulation of an ensemble of electrons can be 
used to assemble a histogram of transit times, which 
in the limit of large numbers of electrons yields the 
transit-time distribution function. This "Monte Carlo" 
method, which incorporates the proper boundary condi- 
tions for our current-based spin detection technique au- 
tomatically, was used to model spin transport through 
doped Si.[7q] 

The technique involves discretely stepping through 
time a duration St, modeling drift with a spatial tr ansla- 
tion v(x)St, and diffusion with a translation of ±\/2D5t 
(where the sign is randomly chosen to model the stochas- 
tic nature of the process), until x > L. The drift velocity 
v(x) can be constructed from a depletion- approximation 
model of the transport layer to determine the electric field 
profile E(x) and a mobility model to get u(-E).[75[ Our 
experimental studies with n-type (« 10 15 cto -3 Phospho- 
rus) 3.3^m-thick Si devices indicate that despite highly 
"non-ohmic" spin transport where the transit time can 
be varied by several orders of magnitude resulting from 
a confining electrostatic potential, Monte-Carlo device 
modeling shows that spin lifetime is not measurably im- 
pacted from its value in intrinsic Si. [76| Fig. [8] shows 
example band diagrams under different bias conditions 
(calculated in the depletion approximation) and corre- 
sponding arrival time distributions from the Monte-Carlo 
simulation. 



VII. DISCUSSION 

There has been significant progress in using ballistic 
hot electron spin injection and detection techniques for 
spin transport studies in Si. However, there are limi- 
tations of these methods. For example, the small bal- 
listic transport transfer ratio is typically no better than 
10~ 3 — 1CT 2 ; the low injection currents and detection 



signals obtained will result in sub-unity gain and limit 
direct applications of these devices. In addition, our re- 
liance on the ability of Schottky barriers to serve as hot 
electron filters presently limits device operation temper- 
atures to approximately 200K (although materials with 
higher Schottky barrier heights could extend this closer 
to room temperature) , and also limits application to only 
non-degenerately doped semiconductors. Carrier freeze- 
out in the n-Si spin detection collector at approximately 
20K introduces a fundamental low-temperature limit as 
well. [§3 

Despite these shortcomings, there are also unique ca- 
pabilities afforded by these methods, such as indepen- 
dent control over internal electric field and injection cur- 
rent, and spectroscopic control over the injection energy 
level [9lj|. Unlike, for instance, optical techniques, other 
semiconductor materials should be equally well suited 
to study with these methods. The purpose is to use 
these devices as tools to understand spin transport prop- 
erties for the design of spintronic devices just as the 
Haynes-Shockley experiment enabled the design of elec- 
tronic minority-carrier devices. 

There is still much physics to be done with ballistic hot 
electron spin injection and detection techniques applied 
to Si. Recently, we have adapted our fabrication tech- 
niques to assemble lateral spin transport devices, where 
in particular very long transit lengths [73j and the ef- 
fects of an electrostatic gate to control the proximity to 
a Si/Si02 interface can be investigated. [92j Hopefully, 
more experimental groups will develop the technology 
necessary to compete in this wide-open field. Theorists, 
too, are eagerly invited to address topics such as whether 
this injection technique fully circumvents the "fundamen- 
tal obstacle^because of a remaining Sharvin-like effective 
resistance. |93| or whether it introduces anomalous spin 
dephasing. [89( 
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